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ABSTBACT 


Various  computer-oriented  numerical  schemes  are 
available  to  evaluate  the  velocity  potentials  for 
inviscid  fluid  flow  past  submerged  objects.  The  use  cf 
a surface  distribution  of  sources  of  uniform  strength 
ever  panels  representing  the  object's  surface  requires 
a relatively  fine  grid  to  obtain  accurate  results.  A 
new  variation  of  this  scheme  using  distributed  sources 
of  linearly  varying  strengths  over  triangular  panels 
appears  to  reduce  the  number  of  panels  needed  for  good 
results,  thereby  reducing  computer  time  and  itorage 
requirements.  In  the  present  study,  this  technique  is 
applied  tc  the  problem  of  determining  fluid  forces  and 
dynamic  response  resulting  from  interaction  between  an 
arbitrarily  shaped  object  and  surf ace . waves. 
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I.  INTRODUCTION 


Ocean  waves  interacting  with  large  marine  structures 
produce  substantial  hydrodynamic  forces.  These  forces  are 
experienced  hy  floating  objects  such  as  ships  or  moored 
platforms  as  well  as  submerged  bodies  near  the  free  surface, 
such  as  submarines  at  periscope  depth  or  bottom  mounted 
caissons  in  shallow  water. 

It  is  desirable  from  an  engineering  standpoint  to  be 
able  to  determine  the  wave  forces  and  resulting  dynamic 
response  for  any  large  ocean  structure  under  design.  Until 
recently,  however,  such  determinations  were  not  possible  in 
the  general  case.  Application  of  the  well-known  Mcrison 
eguaticn  to  calculate  wave  forces  is  limited  to  cases  in 
which  structure  dimensions  are  small  compared  to  wave 
length,  as  in  the  case  of  a small  diameter  pile. 

The  usual  analytical  approach  in  dealing  with  the 
hydrodynamic  interaction  of  large  bodies  with  waves  is  based 
on  a mathematical  modeling  of  the  fluid  flow  using  potential 
flow  theory,  which  assumes  the  fluid  to  be  inviscid  and 
irrotational . While  such  an  assumption  is  unrealistic  for 
steady  flew  past  bluff  bodies  (because  of  boundary  layer 
separation  effects),  experiments  indicate  that  it  yields 
valid  results  for  the  case  of  unsteady  flow  past  an  object 
as  produced  by  surface  waves  provided  wave  amplitude  is 
small  compared  to  the  characteristic  dimensions  of  the 
object . 

Numerical  schemes  based  on  the  use  of  the  Green's 
function  (or  source  distribution)  method  which  utilize 
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high-speed  ccmputer  capabilities  have  had  considerable 
success  in  solving  problems  of  wave/body  interactions. 
Garrison,  et.  al.  [Refs.  1-5]  have  developed  numerical 
procedures  for  calculation  of  wave  interaction  with  large, 
three-dimensicnal,  fixed  and  floating  bodies.  These 
procedures,  which  represent  the  immersed  surface  by  a grid 
of  distributed  source  panels  of  uniform  strengths,  have 
proven  to  be  powerful  in  practical  application  as 
demonstrated  in  Refs.  6-9. 


The  most  severe  limitation  of  the  above  mentioned 
procedures  is  the  rather  large  amount  of  computer  (C.P.U.) 
time  reguired  to  obtain  an  accurare  solution  . Thus,  any 
improvements  to  the  method  which  would  improve  convergence 
of  the  solution  would  be  of  considerable  value. 


A natural  extension  of  the  technique  based  on  a uniform 
distribution  of  sources  over  each  grid  panel  is  to  allow  the 
source  strength  to  vary  linearly  over  the  extent  of  each 
panel.  Since  this  is  not  convenient  when  using  quadrilateral 
panels,  it  is  appropriate  to  divide  the  source  surface  into 
a grid  of  triangular  panels.  The  linearly  varying  source 
strength  may  then  be  represented  in  terms  of  the  values  of 
the  source  strength  at  the  corners  of  the  triangles. 

Webster  [Ref.  10]  has  developed  such  a scheme  for  the 
case  of  uniform  flow  of  an  unbounded  fluid  past  a 
three-dimensional  body.  The  body  was  represented  by  use  of 
the  three-dimensional  source  potential  of  the  form  1/R  where 
R denotes  the  distance  from  the  source.  In  the  present  study 
this  general  method  is  extended  to  the  problem  of 
representing  a body  of  arbitrary  shape  located  in  fluid  of 
finite  depth  in  the  presence  of  a free  surface. 


The  basic  approach  used  in  solving  the  problem  of  wave 
interaction  with  large  floating  objects  is  to  decompcse  it 
into  compcnent  parts  which  are  more  simply  analyzed,  and 
then  superimpcse  the  results  to  yield  the  overall  solution. 
A group  cf  bcundar y- value  problems  are  generated  in  terms  of 
velocity  potential;  one  problem  is  solved  for  each  of  the 
sir  possible  degrees  of  freedom  of  body  motion,  and  one  for 
the  scattering  of  the  incident  wave  by  the  body.  The 
solutions  are  then  formulated  in  terms  cf  the  Green’s 
function.  The  resulting  integral  equation  is  then  converted 
to  a system  of  linear  equations  which  are  subsequently 
solved  by  ccaputer. 


II.  fOS^lATION  OF  THE  HAVE  INTERACTION  PS03I|ji 


A.  fiESTSIClIONS  AND  AS30SPTI0NS 


All  of  the  usual  simplifying  assumptions  incorporated  in 
linear  gravity  wave  theory  are  adopted  in  the  present  study. 
In  summary,  these  are: 

a The  flow  is  considered  to  he  incom pressibie, 
irrotaticnal , and  inviscid. 

a The  wave  amplitude  to  wave  length  ratio  is 

small . 


and  the  dens 

a 

horizontal, 

a 

equation  is 


The  pressure  is  uniform 
ity  is  uniform  throughout 
The  bottom  boundary 
impermeable  plane. 

The  velocity-sguared 
neglected . 

The  waves  are  regular  and 


above  the  free  surface, 
the  fluid. 

is  represented  by  a 
term  in  Bernoulli's 
periodic. 


In  addition,  it  is  assumed  that 

a The  object  has  a smooth, 

surface . 

a The  wave  amplitude  is 
object's  dimensions. 


rigid,  and  impermeable 
small  compared  to  the 


B.  PROBLEM  CESCaiPTION 


The  general  description  of  the  problem  is  shown  in  Fig. 
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1,  in  which  tac  cccrdinate  systems  are  defined.  The  inertial 
coordinate  system,  positioned  at  the  mean  free  surface,  has 
directions  defined  by  x,  y and  z,  while  the  coordinate 
system  positioned  on  the  object  has  directions  indicated  by 
X',  y'  and  z'.  The  origin  of  the  body  coordinates  is  taken 
to  be  at  the  center  of  gravity  of  the  body,  and  the  body  is 
submerged  an  arbitrary  distance,  d,  below  the  mean  free 
surface.  The  incident  wave  is  assumed  to  propagate  in  the 
direction  defined  by  the  angle  JT  indicated  in  Fig.  1. 


The  analysis  of  the  floating  object  located  in  a train 
of  regular  waves  is  decomposed  into  seven  problems  which  are 
uncoupled  and  can  be  solved  independently  and  superimposed 
to  obtain  the  complete  solution. 

The  first  six  problems  are  identical  in  icri,  and 
consist  cf  determining  the  fluid  forces  due  to  oscillation 
of  the  object  in  the  six  degrees  of  freedom,  in  an  otherwise 
still  fluid.  The  six  modes  of  oscillation  are  all  considered 
to  take  place  with  angular  frequency,  0'  = Z7r/T,  which  matches 
the  frequency  cf  the  incident  waves.  These  harmonic  motions 
of  the  object  may  be  expressed  as 


X|,=  Re  [X°  = 


(1) 


o o o 

in  which  X =a 9 when  k=4,  5,  6,  and  where  © denotes  the 
k k k 


angular  displacement  about  the  x',  y'  and  z'  axes. 


respectively.  The  subscripts  k=1,  2,  3 correspond  to  surge, 
heave,  and  sway  modes  of  motion,  respectively,  while  the 
subscripts  k=4,  5,  6 correspond  to  the  rotational  modes  of 
roll,  yaw,  and  pitch,  respectively,  as  depicted  in  Fig.  2. 


o 

The  parameter  X denotes  the  complex  amplitude  of  the  linear 
k 


1 1 


denotes  t.ae  complex 


o 

modes  of  motion,  while  the  term  a 0 

k 

amplitude  of  the  rotational  degrees  of  freedom. 


The  seventh  pronlem  to  be  considered  is  the  interaction 
of  waves  with  the  object  held  fixed  in  position.  This 
interaction  cf  an  incident  wave  with  the  fixed  body  results 
in  "scattered"  waves  which  radiate  from  the  body.  Prom  the 
solution  of  this  problem  the  excitation  forces  and  mcments 
acting  on  the  bojdy  are  determined. 


C.  SQUATICNS  0?  dOTION 


For  a flcatinc  body  with  linear  elastic  constraints  the 
eguation  of  motion  for  the  six  degrees  of  freedom  is  given 
in  Hef.  1 as 
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L-  1,1 


(2) 


where  V - , and  CT  = 21X/t  denotes  the  frequency  of  the 
incident  wave  and  resulting  motion.  (In  Eg.  2 and  ail 
equations  which  follow,  the  use  of  tne  lower  case  letter,  i, 
as  a subscript  denotes  an  index.  Appearance  of  i elsewhere 
in  an  equation  denotes  the  imaginary  number  i=  a/~ 1 ' •)  Tbe 
surface  elevation  of  the  incident  wave  at  the  origin  is 
given  by 


^ = ^°cos  (or-t) 


(3) 


in  which  is  a real  number  denoting  the  half-amplitude  of 
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the  incident  wave.  All  phase  angles  are  referred  to  tne 
wave  crest  with  positive  values  indicating  lag. 


In  Eg.  2,  m denotes  the  dimensionless  body  mass  tensor 

ij 

and  is  defined  by 

^11  “ '^zz  “ "’as  = rn/j>a.^  (4ai 

"’.z=  '"i3  = '"iai^"’zj  = "’3Z.  = ° 

in  which  m denotes  the  mass  of  the  body,  Ji  denotes  the  fluid 
density  and  a denotes  the  characteristic  length  scale  of  the 
body.  Similaily, 

where  the  mciients  of  inertia  are  defined  typically  as 

in  which  the  integration  extends  over  the  complete  mass  of 
tne  body  denoted  by  a. 

The  excitation  force  and  moment  coefficients  are  defined 
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as 


r . = ^.1  s 

J>3  CL^ 


J = t,z^3 


F.  Cmax)  e‘ 


J n V,  5-,  6 


(6) 

(7) 


where  i; denotes  the  phase  shift  of  the  force  with  respect  to 

<J 

the  crest  of  the  incident  wave.  A positive  value  of  Sj 
denotes  a lag.  rhe  term  F (max),  (j=1r  2,  ...6),  denotes 

d 

the  maximum  value  of  the  force  (j=1,  2,  3)  and  moment  (j=4. 


5,  6)  acting  on  the  body. 


In  the  case  of  the  forces  (and  moments)  due  to  the 
motion  of  the  body,  it  is  more  common  to  describe  the 
hydrodynamic  force  in  terms  of  the  added  mass  tensor,  d 

ij 

and  damping  tensor  N . These  parameters  are  defined  as 

ij 


I J 


j^cr^a**Xj 


J J z Jj  2 ^ 


(8) 


(8) 


where  F (max)  denotes  the  amolitude  of  the  force  or  moment 

ij 

in  the  i-th  direction  resulting  from  the  motion  of  the 
floating  tody  in  the  j-th  degree  of  freedom. 

The  forces  defined  in  2gs.  6-9  are  given  in  terms  of 
surface  integrals  of  the  hydrodynamic  pressure  as 


) * J ” I ; 2 J , 4, 


(10) 
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Fi(0  = “ jff  Po7 


J ! ~ ...  (o 


where  F denotes  the  i-th  component  of  wave  excitation  force 
i 

or  moment  and  P denotes  the  i-th  component  of  force  or 

id 

moment  arising  from  the  j-th  component  of  body  morion.  The 


functions  h are  defined  as 
i 


h,  = n. 


^>1=  riy 


hj  = 


-z  n. 

= Z Hx  " X 

= XHy  -Mty)  Hx 


(12a) 


(12b) 


(12c) 


(12d) 


(12e) 


(12f) 


in  which  n.  , n and  n denote  the  three  components  of  the 
X y z 

outward  unit  normal  vector  on  the  surface  of  the  body. 


The  parameters  P , (j=1/  2, 6),  denote  the  pressure 

j 

on  the  immersed  surface  associated  with  motion  in  the  six 

degrees  of  freedom  and  P denotes  the  pressure  associated 

07 


i 


with  the  incident  wave  (0)  and  scattering  (7)  cf  the 


incident  wave.  Finally,  in  Eg.  2 the  dimensionless 

coefficients  K denote  the  hydrostatic  restoration  force 
i ij 

j (or  moment)  tensor  and  are  defined  as 

I 


^zz~ 

s - 

(13a) 

1 z'nyds 

(13b) 

'n^-z'^r\y)ds 

(13c) 

X'z  r\y  ds 

> 

(13d) 

^Zi  ■ 

■ 

f 

1 x'  cis 

(13e) 

[(XV 

is 

(13f) 

where  A denotes  the  waterline  area. 

w 

flooring  line  reactions  depend,  in  general,  on  tne 
mooring  line  configuration,  weight,  shape  (catenary) , 
hydrodynamic  and  elastic  properties.  However,  it  is 
generally  permissible  to  disregard  dynamic  effects  and 
approximate  the  reactions  by  a linear  relationship  with  the 
six  components  of  body  displacements  as  indicated  by  Eg.  2. 
The  term  K*  in  Eg.  2 represents  the  dimensionless  force 

i j 

(or  moment)  produced  in  the  (negative)  i-direction  ny  a unit 

displacement  cf  the  body  in  the  j-direction.  Given  the 
mooring  line  configuration  and  elastic  properties,  the 
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spring  constant  tensor  can  be  evaluated. 


D.  HEPS5.SENIATI0N  BY  VELOCITY  POTENTIALS 


The  fluid  motion  in  each  of  the  seven  flow  situations 
previously  described  can  be  represenred  by  a velocity 
potential.  Specification  of  the  velocity  potential  of  a 
flew  completely  describes  the  flow,  since  fluid  velocities 
at  any  point  are  obtainable  tnrough  the  relation 


where  g(x,y,z,t)  represents  the  velocity  vector, 
represents  the  time  dependent  velocity  potential, 

denotes  the  gradient  of  the  potential.  Thus,  for  the 
six  modes  of  oscillation  of  the  body  in  still  water,  the 
fluid  velccity  vector  is  given  by 

^ = Re.  [v70  (X,y,z)e  (15) 

k=l,z...4 


(14) 


^(x,y,z,t) 


a n 


in  wnich 


4>  (x,y,z) 
k 


denotes  the  complex  space  dependent  part 


of  the  total  potential  induced  by  the  k-th  mode  of 
oscillation.  The  velocity  potential  associated  with  the 
seventh  problem,  waves  interacting  with  the  fixed  body,  is 
actually  composed  of  two  components.  One  component  is  the 
potential  due  to  the  linear  incident  wave,  the  space 
dependent  part  of  wnich  is  given  by 


-if/  f cesh  k(hf-‘/)  \ p ^lhxC0S}frk2.%\n!(]  (16) 

CT  Ccs/1  Kh  I 
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. o 

in  which  “ =K/2  denotes  wave  amplitude  (half-amplitude)  and 
lc=  2TT/l  denotes  wave  number. 

The  second  component  of  the  wave  interaction  potential 
is  due  to  the  presence  of  waves  whose  characterist ics  have 
been  altered  by  interacting  with  the  object.  This  potential 
is  referred  to  as  the  scattering  potential  and  is  denoted  by 
the  subscript  (7) . The  total  potential  for  the  motion  of  the 
body  in  regular  waves  is  expressed  through  the  superposition 


I 


(17) 


E.  THE  ECiJNCARY-V  ALOE  PROBLEM 

Each  cf  the  velocity  potentials  ( S , k=1,2...7)  must 

k 

satisfy  the  continuity  eguation,  which  in  tae  case  of 
irrctaticnal  flow  takes  the  form 


0 , (13) 


In  addition,  the  potentials  must  also  satisfy  tne  kinematic 
bottom  boundary  condition 


|♦K(X,-h,z:)=0 


, K=o,»...7 


(19) 


where  h denotes  the  mean  water  depth,  as  well  as  the 
linearized  dynamic  and  kinematic  free  surface  boundary 
conditions  expressed  by 


(20) 


^k(x,o,z)  _ ^^4>i^(x,o,z)  = o ^ k-Oj\.,.7 

where  g denotes  the  acceleration  of  gravity. 

While  the  free  surface  and  bottom  boundary  conditions 
are  exactly  the  same  for  all  of  the  potentials,  the 
kinematic  boundary  condition  applied  on  the  surface  of  the 
object  is  different  for  each  case.  This  boundary  condition 
results  from  the  assumption  that  the  surface  of  the  object 
is  rigid  and  impermeable;  it  simply  imposes  the  condition 
that  there  be  no  normal  component  of  fluid  velocity  relative 
to  the  surface  at  any  point  on  the  surface,  and  is  expressed 
mathematically  as  follows: 


For  oscillation  of  the  body  in  the  still  fluid. 


-Q  - ~ ^ 


_ Q - 


ny 


- Q - 

Z>n 


C(civy)n2: 


- q = - L<r  e!  Cxn. -(d^y^axl 
'an 


(21a) 

(21b) 

(21c) 

(21d) 

(21e) 

(21f) 


where 

unit 


n , n and  n denote  the  x,  y and  z components  of  the 
X y 2 

outward  normal  vector  on  the  immersed  surface  which  is 
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defined  in  its  aean  position  by  S(x,y,2)=0. 

For  wave  interaction  with  tne  fixed  body, 

COS  tkz  sia  y) 

^ <rco5h(Kh') 

• [r\y  SmKCk  (h^yjI  + i.  cosy +(1^  sin  y)  • Cosh  Ck(hfy)]] 

in  which  y denotes  the  angle  between  the  incident  wave 
propagation  direction  and  the  x'-axis.  Equations  21  and  22 
are  applied  only  on  the  immersed  surface  as  defined  by 
S (x,y,z) =0. 
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III. 


Q.Z  ZM  30UNDA?.Y-VAm;  psoblsm 


A.  GHSEH'S  FONCTIO.’i  SOLUTION 

In  using  Green's  function  representation  of  the  velocity 
potential,  each  of  the  seven  unknown  potentials  is 

considered  to  be  induced  by  fluid  sources  which  are 

distributed  over  the  surface  of  the  body  in  a continuous 
manner.  The  strength  of  the  source  at  any  location  on  the 
body  surface  is  a function  of  position  and  is  denoted  by 
f(^,'^)‘^),  where  indicates  a point  on  the  immersed 

surface.  (Figure  3 illustrates  the  concept  of  the  continuous 
source  d istr ibution . ) The  velocity  potential  associated  with 
each  of  the  boundary  value  problems  is  then  given  by  the 
surface  integral 

(23) 

where  G(x,y,z;  ? ) is  referred  to  as  the  Green's 

function.  Its  value  depends  upon  tne  location  of  the  point 
( > ) on  the  source  surface,  and  the  point  (x,y,z)  in  the 

fluid  field.  The  Green's  function  for  the  present  case, 
which  satisfies  the  bottom  and  free  surface  boundary 

conditions,  is  given  by 

& = */R  + l/R*  + G-*  (2U) 


in  wnich 
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R=V(s-x')^ 


(25a) 


R*  ♦-  (y  taWt^f-K^-zV 


(25b) 


The  G*  term  is  given  by  Hef.  3 as 

^*Cx,y,z;?,7,^)  = Z RV.  f 

''o 

6^  cos/?  Iji  (y^h)]  Jr.  Uir)JjA 

^ smh  (jfh)  - rcosh(^^) 


■h  L ZT^  (k^  -V^)  C0%h[k(7;th)]  cosh  [kCv^h)]Jr.(kr) 

k'^k  -v"^h  -hV 


in  which 


r^Ax-%f  +(z-<r 


The  term  J denotes  the  Bessel  function  of  the  first  Kind, 
0 

cf  order  zero,  and  V represents  dimensionless  wave 

freguency.  E.V.  indicares  the  principal  value  of  rhe 
infinite  integral. 


An  alternate  series  form  of  the  complete  Green's 
function  is  given  ty 


0=  2.TT  [iC  (hr^)]  CosKLk('(tk)JL'(,(Kr)  ~ i 

k^h  -v^h  tV 
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2 ccs^^(„h)]ccs[^^  r) 

lc~  I k 


(28) 


in  which  1 denotes  the  Bessel  function  of  the  second  Icind, 
0 

of  order  zero,  K denotes  the  modified  Bessel  function  of 

0 

the  second  kind,  of  order  zero,  and  represents  the  real 
positive  roots  of  the  equation 


(29) 


In  numerical  evaluations  either  form  of  the  Green's 
function  nay  be  used  with  the  exception  of  the  case  where  r 
approaches  zero.  Here  the  series  form  given  in  Eg.  28  is 
singular  so  the  integral  fora  given  by  Eq.  26,  which  has  a 
well-defined  1/R  singularity,  is  used-  However,  when  r is 
not  small  the  series  form  is  generally  used  for  numerical 
calculation  because  it  requires  less  computer  time  for 
evaluation. 

The  Green's  function  representation  of  the  potential 
given  by  Eg.  24  satisfies  the  laplaca  equation  as  well  as 
the  free  surface  and  bottom  boundary  conditons  for  all 
cases.  The  source  strength  function  f { f is  unknown 

and  is  determined  by  application  of  the  kinematic  boundary 

conditions  on  the  body  surface  (Egs.  21  and  22).  This 
results  in  the  following  integral  equation: 

(x,y,z)  (30) 

O *v 

k=  I,:i  ...7 

in  which  “3  Van  denotes  the  deriviative  of  the  Green's 
function  in  the  outward  normal  direction,  which  is  given  by 


i 


1 

i 

' 

i 


23 


where  q represents  the  outward  unit  normal  vector  at  (x,y,z) 
and  where  (x,y,z)  now  represents  a control  point  on 
S {x,y,z)  =0. 


For  the  seven  problems,  all  quantities  in  Sg.  30  are 
known  except  for  the  source  strength  functions  f ( ^ 9 ), 

(k=1,2,  ...7).  However,  tne  integral  aquation  cannot  be 

solved  analytically  for  an  aroitrarily  shaped  body,  and  a 

numerical  approach  must  be  adopted.  The  numerical  approach 

developed  herein  involves  discretizing  both  the  source 

strength  function  and  the  body  surface  boundary  condition 

functions,  g , (k=  1 , 2 , . . . 7)  . 
k 


3.  DISCRETIZATION  OF  IH2  PR03L2M 


I • Uniform  Source  Strength  Element  Method 


Ihe  technique  previously  used  [Ref.  1]  to  discretize 
the  integral  equations  given  by  Eg.  30  is  to  divide  the 
source  distribution  into  a grid  of  N quadrilateral  panels, 
and  to  assume  that  the  source  strength  has  a distinct, 
constant  value  over  each  panel.  Hence,  for  any  of  the  seven 

flew  problems  of  interest  here,  f denotes  the  strength  of 

''j 

the  source  everywhere  on  panel  j,  (j=1,2...N),  for  the  k-th 
problem  (k=1,2  ...7).  The  area  of  panel  j is  denoted  ty 
ASj  and  the  induced  potential  anywhere  in  the  fluid  thus 
results  frem  a summation  of  the  effects  of  all  the  source 
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panels.  Figure  4 depicts  the  discretized  source  surface  and 
resulting  potential. 


Additionally,  the  noraal  boundary  condition 

functions  g (x,y,z),  (k=1,2...7),  of  Eg.  30  are  discretized 

k 

ny  applying  the  kinematic  body  surface  boundary  condition  at 


only  M selected  points  on  the  surface.  In  practice,  the 
number  of  boundary  condition  ("control")  points  are  chosen 
to  be  egual  (M=N) . For  convenience,  the  control  points  are 
chosen  to  lie  at  the  centroids  of  the  source  panels.  Thus, 
for  each  discretized  integral  eguation,  a system  of  11  linear 
equations  in  N unknowns  results: 


ziSj 


i,Z...7 

L-I,Z...N  (31) 


in  which  (x  ,y  ,z  ) represents  a control  point  i on  the  bodv 
i i i 

surface,  at  which  the  boundary  condition  is  being  applied, 
and  g^  denotes  the  value  of  g^  (x,y,z)  at  point  i.  Eguation 


31  is  expressed  in  matrix  form  as 


W " 1,2...  7 (32) 

where  [o<j  is  an  !<xN  matrix  composed  of  elements 


strengths  and  the  kinematic  noundary  condition  applied  at 
the  N centred  points,  respectively,  for  the  k-th  flow 
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problem.  It  should  be  noted  that  the  matrix 
only  on  the  body  configuration,  water  depth, 
the  motion  ard  is  the  same  for  all  seven  cases 


[ o<  ] depends 
and  period  of 


Ihe  discretized  potential  which  is  induced  by  the  N 
source  panels  is  specified  by  a matrix  equation  similar  to 
Eg.  20: 


{%]  = [^j  1^.] 


The  martix  [f3]  is  composed  of  elements 


(35) 


As  in  the  case  of  matrix  [o<],  matrix  [/3]  is  the  same  fcr  all 
seven  flew  cases. 

Integration  of  the  Green's  function  and  its 
derivatives  in  Eqs.  33  and  35  poses  no  great  difficulty.  It 
is  most  convenient  to  separately  integrate  the  singular 
terms.  Special  ' considerations  in  performing  these 
integrations  and  necessary  formulas  are  provided  in  Ref.  2. 


A computer  solution  of  Eg.  32  is  accomplisned  using 
matrix  inversicn  and  multiplication,  yielding  the  previously 
unknown  values  of  the  source  strengths.  Values  of  velocity 
potential  at  the  control  points  (panel  centroids)  are  then 
obtained  directly  using  Eqs.  34  and  35.  Eluid  velocities  at 
the  control  points  are  similarly  obtained  after  first 
differentiating  the  elements  of  the  [jQ  ] matrix  with  respect 
to  the  x,y  and  z directions.  Dynamic  pressures  are 
calculated  from  the  linearized  form  of  Bernoulli's  equation. 
Finally,  the  forces  acting  on  the  object  are  determined  by 
pressure  integraticn  (or  more  precisely,  summation)  over  the 
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surface  area  of  the  object  as  indicated  by  Eqs.  10  and  11. 


1 

t 

1 

2*  Method  of  Linearly  Varying  Source  Strength  Elements  i 

i 

An  alternative  method  of  discretizing  tne  Green's 
function  integral  equations,  proposed  by  Webster  [3ef.  10] 

for  uniform  unbounded  flow  past  a submerged  body,  utilizes  j 

triangular  source  panels  over  which  the  source  strength 
function  is  considered  to  vary  linearly.  In  this  way,  the  ; 

discontinuities  or  "jumps"  in  source  strengths  at  adjoining 
panel  edges  (which  existed  in  the  previous  scheme)  are 
eliminated.  Also,  since  source  strengths  are  not  held 
constant  over  the  extent  of  each  panel,  the  discretization 
results  in  a less  crude  approximation  of  the  continuous 
strength  funotion.  An  additional  "smoothing"  effect  on  the 
velocities  induced  near  tae  body  is  achieved  by  submerging 
the  source  panel  surface  inward  a selected  distance  from  the 
actual  body  surface,  although  the  kinematic  body  boundary 
condition  is  still  applied  at  points  on  the  true  body 
surface.  It  is  noted  that  while  some  smoothing  effect  can 
be  accomplished  by  submerging  the  sources,  this  is  not  a 
necessity  as  long  as  the  singularity  is  accounted  for 
properly . 


rigure  5 

dap 

icts 

the 

linearly 

varying 

source 

strength  function 

for 

a single 

triangular 

panel. 

The 

local 

coordinate  system 

for 

each  panel 

is  defined 

such 

tha 

t two 

triangle  corners 

(de 

noted 

by 

a and  d) 

lie  on 

the 

local 

x-axis,  and  the  t 

bird 

corner 

(b)  lies 

on  the 

po 

sitive 

y-axis.  Appendix  A gives  the  transformation  between  the 
body  coordinate  system  and  the  local  coordinate  system  for  a 
given  panel. 

In  panel  coordinates,  the  linearly  varying  source 
strength  function  may  be  written 
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(36) 


¥ 


in  which  ®<  ^ /3  and  'If  denote  constant  coefficients  and  the 
tildes  indicate  local  coordinate  values  as  defined  in  fig. 
5.  The  values  of  the  source  strength  at  the  three  corners 


are 

1 

1 

III  ’ 

1 

? i 

, 

1 

{ 

I 

I 

i 

1 

( 

1 

(37a) 

'3 

o' 

III 

(37b) 

h = fca.o) 

(37c) 

constant 
triangle 
corners : 

Substitutuion  of  Egs-  37  into  Eg.  36  allows  the 
ooefficients  to  be  expressed  in  terms  cf  the 
geometry  and  the  values  of  the  strength  at  the 

(38a) 

(38b) 

- [(a-a)ft  tJ 

(38c) 

Thus,  the  contribution  to  the  potential  at  (x,y,z)  due  to 
the  single  panel  of  area  is  denoted  by  A^(^/y»z)  and  is 
given  by: 
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ff  c/s 


(39) 


wh<=re  ^ and  denote  local  panel  coordinates  and  (x,y,z) 
denotes  a general  point  expressed  in  terms  of  local 
coordinates. 


Substituting  Egs.  38  into  Eg.  39  and  collecting 
terms  alleys  the  contribution  to  rhe  perential  to  be 
expressed  in  terms  of  the  values  of  the  source  strength  at 
the  cornets  of  the  panel: 


in  which 

A</>t  = J G-CCf ; ^ ?)  <^s 
A </>J  =1 

''as 

For  convenience  and  will  be  referred  to  as 

simply  the  "corner  integrals"  associated  with  corners  a,b 
and  d,  respectively,  for  a given  triangular  panel. 

Ihe  velocity  induced  by  a single  source  panel  is 
obtained  by  taking  the  gradient  of  Eg.  aO: 


(40) 

(41a) 

(41b) 

(41c) 


A ^ ( x~,  Y,  z)  (A  4)^)  ^ f V(A  dij) 
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The  gradients  of  the  corner  integrals  are  given  by 


^ = jj  ^ ^ ^ 

= I..  TsSf 


(43a) 

(43b) 

(43c) 


The  evaluation  of  the  three  corner  integrals  (Eg. 
41)  and  their  gradients  (Eg.  43)  must  be  accomplished  in 
local  coordinates  for  each  panel.  In  perrorainc  the 
integrations,  it  is  most  efficient  to  separate  the  Green's 
function  into  its  singular  and  non-singular  components,  and 
treat  the  integrals  separately,  as  follows: 

- II  Ki  ( <-  \{  K ^ j( 

hiS  ''■’a.s  ^•'aS 

' fLi  ^ t j[  ^4  * jjj'i 

WdC/R.Vsfjj'  Wi('/R')ds 

V(A(ti„Vf  W,7(I/r)«Is +((  W^Vfl/R')ds  + ffw<^(7(&*)<yi 

'JaS  '-'oS 


(44a) 

(44b) 

(44C) 
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in  which  ii 


# 


K , and  W are  linear  functions  in  ? and ^ and 
b d ^ 


are  given  by 


Wa'-- 

bd  t b f 

*f7r 

(46a) 

-n 

(46b) 

UK  b 

II 

bfl  — 'o  ? - a ^ 

(46c) 

Affendix  B provides  foraulas  developed  by  Webster 
[Hef.  10]  fcr  the  closed  fora  evaluation  of  those  integrals 
in  Eqs.  44  and  45  which  involve  derivatives  of  1/R  and  1/R' 
as  well  as  siailar  foraulas  developed  by  Yeung  and  Bai  [Ref. 
11]  for  integrals  of  1/R  and  1/R*.  The  integrands  are 
singular,  and  it  is  therefore  necessary  to  evaluate  the 
integrals  analytically. 

Integration  of  the  G*  components  of  Egs.  44  and  45 
is  made  possible  by  assuming  that  the  function  and  its 
gradient  vary  linearly  over  the  triangular  panel,  (as  was 
done  for  the  source  strength  function) . This  is  a 
reasonable  approximation  since  G*  is  a well-behaved  function 
that  varies  slowly  with  position  on  the  body  surface.  The 
wave  length  of  its  oscillation  is  approximately  the  same  as 
the  wave  length  of  the  incident  wave,  and  this  tends  to  be 
large  relative  to  the  panel  dinensions  in  typical 
applications . 

With  G*  and  VG*  expressed  as  linear  functions  of 
their  values  at  the  panel  corners,  the  terms  in  Egs.  44  and 
45  involving  these  two  functions  may  be  expressed  as 
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(47a) 


- I-a  (S-* 

(47a) 

= rja,b,<i, 

(47b) 

Wj  s = Zii a, b, J,  &y , Q* , g/) 

(47c) 

(48a) 

Wb  v&'ds  = rj  U,h,I,  xjQfyG.l) 

(48b) 

(48c) 

where  I , I , I , I ' , I ' , I ' , represent  the  integrals  cf  the 
a b d a b d 

product  cf  the  two  linear  functions  over  the  triangular 
area.  The  integrals  in  Egs.  47  and  48  were  evaluated  using 
the  formulas  given  in  appendix  C for  integrating  the  product 
cf  two  linear  functions  over  a triangular  area  and  the 
resulting  integrals  are  given  in  Appendix  D 

The  total  potential  and  velocity  at  any  point  in  tue 
fluid  is  a result  of  the  combined  effects  cf  all  the  source 
panels,  for  a grid  composed  of  a panels,  the  potential  is 
given  by 


(49) 


and  the  velocity  vector  is  given  by 


q(x,,,i)=E  '^"1 

X k k k k k 
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in 

which 

, f , denote  the 

strengths  at 

the  three 

corners  of 

the  k-tn 

panel . 

It 

should 

be  emphasized  that  the  corner 

integrals 

and 

their 

gradients 

are  evaluated  in 

their  own 

indi widual 

panel  coordinate  systeas.  Iransf oraation  back  zo  the  global 
coordinates,  described  in  Appendix  A,  aust  be  accoaplished 
prior  to  performing  the  suafflations  of  Eqs.  49  and  50. 

Ihe  sumaations  given  in  Eqs.  49  and  50  can  be 
re-arranged  by  reversing  the  order  of  rhe  suaaation.  That 
is,  if  coefficients  of  the  strength  at  a given  node  are 
first  suaaed,  then  the  suamation  aay  ne  written 


= E Tj  A(t); 
j=i 

(51) 

E fj  9'(A4>,(x,y,Z)) 

(52) 

in  which  f denotes  the  source  strength  at  corner  node  j 
3 

(j=1,2...N)  and  A represents  the  sun  of  the  corner 

3 

integrals  associated  with  node  j,  for  all  panels  which  share 

node  j as  a mutual  corner.  A grapnic  representation  of  the 
meaning  of  AcJ)  , is  provided  in  Fig.  6 for  clarity,  since  it 

is  difficult  to  express  mathematically.  As  indicated  in 

Fig.  6 A(i)  denotes  the  contribution  to  the 

3 

potential  at  (x,y,z)  from  the  panels  surrounding  node  j. 

Each  of  the  surrounding  panels  may  be  considered  to  have 
unit  strength  at  the  corner  coincident  with  node  j and  zero 
at  the  other  two  corners  as  depicted  in  Fig.  6. 
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The  kinematic  boundary  condition  at  the  body  surface 
may  now  he  applied  at  N selected  points.  It  is  convenient 
to  use  points  on  the  surface  which  lie  directly  outward  from 
the  corner  nodes  of  the  submerged  source  surface.  In 

t 

practice,  it  is  best  to  initially  define  the  nodes  on  the 
actual  body  surface,  and  then  establish  the  points  on  the 
submerged  grid  at  seme  pre-selected  submergence  distance 
along  the  corresponding  inward  normal  vectors. 

Application  of  the  kinematic  boundary  condirions  for 
each  of  the  seven  problems  of  interest  results  in  the 
discretized  version  of  the  integral  equation: 


1)2.. .-.Ij  1-1^2... .M 


where  the  repeared  index  j indicates  summation.  The  term 
denotes  the  sum  of  the  corner  integrals  for  node  j, 
evaluated  at  control  point  i.  The  derivative  with  respect 
to  the  outward  normal  vector  is  obtained  by 


where  n is  the  unit  outward  normal  vector  at  tae  control 
i 

point  i,  on  the  body  surface. 


Equation  53  may  be  expressed  more  simply  as 


t-j  ly  • i\! 


, l<  - .. . 7 


where  the  NxN  matrix,  [oi  ],  is  composed  of  elements 


(56) 


Equation  55  is  solved  by  computer,  (for  k.  = 1,2...7) 

using  matrix  inversion  and  multiplication,  yielding  the 

unknown  source  strengths.  Substitution  of  the  source 

strengths  into  Egs.  51  and  52  with  evaluation  at  (x  ,y  ,z  ) 

i i i 

(i  = 1 , 2 , . . . N)  , gives  the  values  of  the  potential  and  the 
fluid  velocities  at  the  control  points  on  the  body  surface. 

Ihe  dynamic  fluid  pressure  is  determined  using  the 
linearized  form  of  Bernoulli's  equation.  For  oscillation  of 
the  body  in  its  six  degrees  of  freedom,  the  pressures  at  the 
surface  nodes  (denoted  by  i)  are  given  by: 


II 

Re 

» — 1 
-p 

b 

...7 

(57) 

Po^• 

= Re  LiPcr  (4>o 

4- 

(58) 

where  jO 

denotes  fluid 

density. 

and  P 

and 

denote 

dynamic 

pressure  and 

potential 

, respect 

ively. 

at 

node  i 

(denoted 

by  the  subscrip 

t i)  due 

to  motion 

in  the 

k-th 

mode 

of  oscillation. 


The  pressure  integrals  given  by  Eqs.  10  and  11 
cannot  be  utilized  directly  to  obtain  total  forces  and 
moments,  because  only  discrete  values  of  pressures  are 
available.  However,  it  is  consistent  with  the  numerical 
scheme  to  assume  that  rhe  pressure  varies  linearly  over  each 
triangular  panel.  Based  on  this  assumption,  the  pressure  at 
any  point  on  a surface  panel  is  expressed  as  a linear 
function  of  the  panel  geometry  and  the  values  at  the  corner 
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I 


nodes,  which  are  given  by  Egs.  57  and  53. 

The  forces  and  moments  are  then  calculated  using 

^ E j(  Pk  hi  <”l 

[[  Poyh/c/s  , 

where  F denotes  the  i-th  component  of  wave  excitation  force 
i 

or  moment  and  F denotes  the  i-th  component  of  force  or 

ik 

moment  arising  from  the  k-th  component  of  body  motion.  The 

summations  of  Eqs.  59  and  60  are  carried  out  over  the  M 

panels  on  the  actual  body  surface.  The  functions  h are 

i 

given  by  Eg.  12  and  are  functions  of  position  on  the 

immersed  surface.  It  appears  to  be  consistent  with  the 
numerical  scheme  to  approximate  these  functions  as  linear 
functions  so  that  they  may  be  defined  at  interior  points  by 
their  values  at  the  corners  of  the  triangles.  Thus  the 
integrands  of  Eqs.  59  and  60  contain  the  products  of  two 
linear  functions.  The  integrals  are  evaluated  in  a 
straight-forward  manner  using  the  integration  formulas  given 
in  Appendix  C . 

The  resulting  forces  F (i ,k  = 1 , 2 . . . 6 ) and  F 

ik  i 

(i=1,2...6)  are  substituted  into  Eqs.  6,  7,  3 and  9 to 
evaluate  the  excitation  force  coefficients  C ( j= 1 , 2, . . . 6)  , 

j 

and  the  added  mass  and  damping  coefficients,  K , and  N , 

ij 

(i, j=  1, 2, . . . 6)  . Finally,  the  equations  of  motion  for  the  six 
degrees  of  freedom  (Eq.  2)  are  applied  to  determine  the 
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IV. 


N'JtMSRICAL  RE^LTS 


A conputsr  code  based  on  the  nuaerical  method  outlined 
was  developed  to  calculate  nhe  hydrodynamic  coefficients  as 
well  as  the  dynamic  response  of  a floating  body  in  waves. 
Using  this  program,  example  calculations  were  made  for  two 
simple  geometric  shapes;  a vertical  circular  cylinder  which 
extends  from  the  bottom  and  passes  through  the  free  surface, 
and  a semi-immersed  sphere.  Computer  produced  drawings  of 
the  grid  configurations  used  are  shown  in  Pigs.  7-10. 

The  vertical  circular  cylinder  represents  the  only 
three-dimensional  geometry  for  which  a closed-fora  solution 
exists  (of  the  type  of  interest  in  the  present  study),  and 
therefore,  it  is  of  interest  for  making  a comparison  of 
numerical  results  computed  by  the  present  method.  HacCamy 
and  Fuchs  [Ref.  12]  have  developed  a closed  form  solution 
for  the  horizontal  force  acting  on  the  cylinder  and  this 
result  is  plotted  in  in  Fig.  11  in  the  fora  of  the 
dimensionless  amplitude  of  the  force  versus  the  wave  length 
parameter,  27ra/L, 

Corresponding  results  computed  by  the  triangular  panel 
method  using  grids  of  two  different  finenesses  are  shown  on 
Fig.  11  for  comparison  with  the  results  of  MacCaay  and 
Fuchs.  These  results  indicate  rather  rapid  convergence  of 
the  solution,  particularily  at  small  values  of  2TTa/L.  The 
slower  ccnvercence  (or  greater  error)  at  the  larger  values 
of  2 n a/L  apparently  indicates  the  inaccuracies  associated 
with  the  assumption  of  linearly  varying  source  strength.  As 
2 7T  a/L  becomes  large  (the  wave  length  becomes  small)  tne 
variation  of  the  G*  part  of  the  Green's  function  with 
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A 


\ 

\ 

\ 

\ 

; distance  beccmes  more  rapid.  Thus,  as  2TTa/L  increases,  a 

' point  is  reached  where  the  assumption  of  the  linear 

' variation  ci  the  Green's  function  becomes  apprecianly 

invalid.  Of  course,  as  the  fineness  of  the  grid  is 
i increased  the  numerical  results  tend  to  agree  witn  the 

j closed-form  exact  results  up  to  higher  values  of  2i1a/L. 

; Havelock  [Hef.  13]  has  presented  results  for  the  added 

mass  and  damping  coefficient  in  heave  for  a semi-immersed 

sphere  and  these  results  are  presented  in  Fig.  12  for 

I comparison  with  numerical  results  computed  by  the  present 

t method.  The  figure  snows  the  heave  added  mass  and  damping 

coefficients,  M and  '3  , respectively,  for  the 

22  22 

semi-immersed  sphere  along  with  the  numerical  results 

computed  using  grids  of  two  different  finenesses.  These 
results  indicate  that  the  numerical  results  are  tending  to 
converge  although  even  for  the  finer  grid  the  results  are 
not  yet  converged. 

Moreover,  given  the  damping  coefficient,  the  wave  force 
coefficient  can  be  computed  by  use  of  the  Haskind's 
relations  which,  in  general,  relate  the  force  coefficients 
to  the  iar-field  solutions  of  the  corr esponding  radiation 
problem.  For  the  case  of  the  heaving  motion  or  an 
a xi-syam etric  body  a closed  form  relationship  exists  between 
the  heave  damping  coefficient  and  heave  excitation  force. 
Thus,  Havelock's  results  for  the  damping  coefficient  in 
heave  were  used  to  compute  the  neave  excitation  force  and 
this  result  is  presented  in  Fig.  13  for  comparison  with  the 
numerical  results. 

The  numerical  results  presented  for  comparison  with 
Havelock's  results  in  Fig.  13  show  a trend  toward 
convergence  and,  in  fact,  the  results  corresonding  to  the 
finer  grid  indicate  adequate  agreement. 
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Finally,  the  dynamic  response  in  heave  and  surge  for  the 
case  of  the  semi-iaaersed  sphere  are  shown  in  Fig.  14. 
These  results,  showing  the  dimensionless  response  in  terms 
of  the  response  amplitude  to  wave  amplitude  ratio,  indicate 
very  little  difference  between  the  the  two  grid  finenesses. 


0 


V. 


CONCLUSIONS 


A new  nuierical  procedure  based  on  the  use  of  triangular 
source  panels  has  been  developed  for  computing  the 
interaction  of  fixed  or  floating  bodies  with  waves. 
Numerical  results  based  on  this  procedure  compare  well  with 
existing  results,  and,  in  general,  the  procedure  appears  to 
converge  to  the  correct  solution. 

The  value  of  the  triangular  panel  procedure  as  coapared 
to  the  procedure  based  on  quadrilateral  panels  of  uniform 
strength  remains  uncertain.  This  will  be  determined  only 
through  experience  in  application  of  the  two  methods. 
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APPENDIX  A 


COORDINATE  TRANSFORMATION 


In  the  following,  the  coordinate  transfor nation  is 
developed  which  relates  the  local  coordinates  atrached  ro  a 
particular  panel,  as  indicated  in  figure  la,  to  the  gichal 
coordinates  with  origin  located  at  the  mean  water  level.  The 
global  system  is  denoted  by  0(x,y,z)  and  the  local 
coordinate  system  is  denoted  by  0(x,y,z). 


Figure  la-Coordinate  System  Locations 


The  location  of  the  origin  of  the  local  coordinate 
system  expressed  in  global  coordinates  is  given  by 

[Xo,yo,Zc]  = [(*d  vxj)),  {rA 


the  glchal 


in  which  (x  ,y  ,z  ) , (x  ,y  ,z  ) and  (x  ,y  /Z  ) are 
aaa  bbb  ddd 

system  cocrdinates  of  the  triangle  corners  a,b  and  d, 

respectively,  and  where 

o<  — ( ^ i'^a. 

V iVa-Yd)^ 

The  unit  vectors  in  the  x,y  and  z directions  of  the 
local  coordinate  system  are  given  in  terms  cf  their 
respective  global  system  components  as 

i ^ ^^JTj  A 

^ ^ k 

l ^ ^ 

in  which  i,  j and  k are  the  unit  vectors  in  the  x,  y and  z 
directions,  and  i,  J and  k are  the  unit  vectors  in  the  x,  y 
and  z directions. 

The  terms  u , u , ...u  are  defined  as  fellows: 

'i  *1  s 
^K,  = 

Ux^=  Tal 

U/3= 

Uy,  = (xi-X,)/l^fel 
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1^2.3=  LViU,,-  u,,ay_j 
Uz3  = U,.U,_^-U<3U,, 

where 

1 cia\  t ' 

and 


l°bl  =y(X,-XO^  ^(Y»-?^,)'‘  fCz^-Zb)'' 


The  transformation  matrix  is  then  defined  as 


Transformation  between  global  coordinates  ana  local 
coordinates  is  accomplished  in  the  following  manner: 


fx' 

’x-x,l  1 

< y 

X 

h 

z' 

V. 

= [T]  . 

y - Xc  • 

^ > 

l = [Tf  y 

Z 

> 

+ 

Vo 

Zc 

1 / 
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J 


in  which  [T]  denotes  the  transpose  of  the  aiatrix  [T]. 


Similarly,  a vector  quantity,  F,  evaluated  in  the  local 
coordinate  system  is  transformed  bacit  to  global  coordinates 
using 


in  which  T , ? and  F are  the  x,  y and  z components  in 
X y z 

local  coordinates  and  F , F and  F are  the  ccrrespcnaing 

X y z 

components  in  the  global  system. 


APPENDIX  B 


FOPaULAS  FOH  EVALUATING  INTEGRALS  INVOLVING  1/R  V2’ 


The  terms  or  Eg.  44  which  involve  1/R  and  1/5'  are 
evaluated  as  follows,  using  the  panel  geometry  defined  by 
Fig.  It. 


Figure  It-  Panel  Geometry 


(V3^t-ad)Fab/c^  - Fjt,  ^ [w’-Cb"''-‘iclK2^h^"')]l^^/c 


+ [2y/-d(y^-e5^)]loS 


p 


Jba-Va  ^5  __  _j r _ A 

^ ^nb(a-cl)R  zb{a-d)ii‘n  f- ^ 

+(b^+ad)Fjb/<^^  - f (b‘^ad)(z^f 
^ Uyp-a(yS2")J  X^a  ] 
ll^^rTbR*^^  - T^xx[2.yz  f + R^  ^aF^ /c^ 

in  which  the  following  definitions  are  applicable 


rl-V(X-AX  r y’-  r£"' 
(^-  b)  ^ ^ Z^' 

-V6r-J)"+  r 


/5=:  bX  -ray  -ah 

A ^ -^r 

= - by  tly'^ 


a ^/FTtF 

a ~y 0-  ^ ^ b^' 

- h%  f ch  y - Jb 
--ci^  -hr,  -y* 

r^x  - bj  f-b^ 


A 


+ r„„-' [-l!ZzlA]  - Un" 

l<\  f^l  ^ '-y  J_  — ' 


Zb  r. 


-^?.7  - ^oa  Lrx  z (x-a)]/[r^  zCx-d'}'] 
^ ^ i A)1 

//,  - e ra  Z.f’.Vre  r,  r//)] 


! 


47 


48 


((  b^-b?-£i  ^ --  ( - . i ^ ') 

j)^^  HiTb(a-<-i)  ^ ^ i b(a-d) 


T d^‘>!  _a£l—^^rr  - ^ ) 

Ls L a 3 tyo(cl-J )*^  bC(^<^-:OJ 


(X  -ci 


£3  hei^-'-T)' 

It.  If  (^-0  b 


r hP  - [H  ^ 

^ L C*  b c.  .]  -^Uk>  ^ 


o <=- 


; i^it] 


— hd  t ( ^/q\S  S 

HVbQi-d) 


b(a 


~dzl/j  ] 


(L^S&Tf  y '*>*  = y ^-<^y;'^^; 

+cziii 

j{  4{y  T'f 

AS 
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APPENDIX  C 


PCHMULAS  FOR  INTEGRATION  OF  LINEAR  FUNCTIONS  OVER  TRIANGULAR 

AREAS 


All  rormulad  given  are  appiicaale  for  a rriangle 
positioned  as  shown  in  Fig.  1c,  with  a>d  and  b>0. 


Figure  1c-  Triangle  Geometry 

If  f(x,y)  and  g(x,y)  are  any  linear  functions  expressed 
as 

f(x.y)  = o<  + (3  X 1-  ^ 7 
g(X,7)=  o<'  r/3'x 

and  the  function  values  at  the  triangle  corners  are  denoted 


’I 


by  £ , f , f , g , g , g , then; 
a h d a b d 


of  r ^ f'a. 

a -d 


(3- 

a~d 

X ~(<^~d)  rdf^ 

bid-ti) 


3a 

b(d-t) 


and  the  integrals  of  f(X/y)ds  and  f (x,y)  g (x,y)  ds  are  given 
by 


fi  f(x,/)ds  ^ w, 

S 3 

jj  f C^>/)  0 ( ^)  y)^-^  = o<cx' w,  t-(pc<'  to<3')  W; 


and 


+ (<X  X ¥o(  XOW3+  (Y/GV  G <y*)  Wl(.  t /3((9  ' W_f  +"  ^ 


where 


w,= 

\\ 

Wj=  ' 
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APPENJIX  D 


INTSGHATICN  FORMULAS  FOR  INTEGRALS  INVOLVING  G*  ANL^VG* 


The  assumption  of  a linear  variation  of  G*  ( ; x, y , z) 
over  a given  triangular  panel,  and  the  suoseguent 
application  of  the  formula  for  the  integral  of  the  product 
of  two  linear  fu notions,  given  in  Appendix  C,  leads  to  the 
following  expressions  for  Sgs.  47  and  46  which  involve  the 
G*  term.  (refer  to  Fig.  1o  for  applicable  triangle 
gecmetry . ) 


ff  -bd  fb 

))  hnb{d-'j) 


G'^ds  - 


+ C ^ A ^ -t-  A ^1 


A 77  [a~d) 


r7tr[^^,E, 


B 


G*  ds 


/ 

MTl  h(j-J) 


1-c 


in  which 

17,  ^ G-’  - 

‘ a <d 

U&-; 
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A,-- 

b (za-d) 

A.= 

:z^ 

B,  - 

(a  b 

6.= 

i z 

B3  = % 

D,= 

b (zd  ta) 

Pi  - 

O3  ■=  % 

Xlf 

ZH 

Ihe 

integratio  ns 

involving 

the 

gradients  of  G-  similarly 

result  in  the  following  expressions: 


1 fCVC^)A^f(vQAj 

'aS  ^ 


rr  bgrh  i-ci'i 
\s 


t 

h'il  b 


V&*ds 


^ 7f  b(ii  -J) 


where 


— I ^Jk.  ^ j ^ j,-  i)Cr  (a  ,<^) 

‘ ^x  -ay  3Z 

— I ^C^^CdjO)  _ y ^ 

SX  ^ dz 

vc^=  [(a-.n(r 
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n{  i J f-  k 

\ 7\x  ■?!'.  'hZ  / 


. j /b 

^ ;)x  ^-2:  ^-1 


y c — {jl  ( L y.  i ^’V^o')  1^  '3«p  %)0)\ 

^ 2X  ^ / 


- rr  irl''’-’’'  ^ i i£;6/)  ^ ^v^.. 

-3= 


Pig.-l-  Problem  Definition 
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Fig.  6-  Contribution  to  Potential  at  Node  i 
From  Panels  Surrounding  Node  z 


L 
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Figure  10  - 


1 


Pig.  13-  Heave  Excitation  Porce  for  a Ploatine- 
Hemisphere 
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